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[ij  This  work  discusses  potential  transport  pathways  of 
terrigenous  material  in  the  Gulf  of  Papua  (GoP),  New 
Guinea,  using  Lagrangian  tracers  as  proxies  for  clay 
minerals  and  finer  particles.  The  tracers  are  transported  by 
currents  from  the  East  Asia  Seas  implementation  of  the 
Navy  Coastal  Ocean  Model.  The  results  suggest  that  clay 
minerals  input  from  rivers  along  the  northwest  coast  of  the 
GoP  accumulate  on  the  inner  shelf  with  high  concentrations 
near  the  mouth  of  the  Fly  and  Purari  Rivers.  Finer  particles, 
which  are  also  representative  of  dissolved  metals,  arc 
transported  eastward  along  the  GoP  coast  into  the  Solomon 
Sea  as  well  as  into  the  Torres  Strait  to  the  west.  The  results 
also  suggest  that  some  finer  particles  are  entrained  by  eddies 
within  the  northern  Coral  Sea.  These  results  are  in 
qualitative  agreement  with  observations  from  the  region. 
Citation:  Keen,  T.  R.,  D.  S.  Ko,  R.  L.  Slingerland,  S.  Riedlinger, 
and  P.  Flynn  (2006),  Potential  transport  pathways  of  terrigenous 
material  in  the  Gulf  of  Papua,  Geophys.  Res.  Lett .,  JJ,  L04608, 
doi:  1 0. 1 029/2005GL0254 1 6. 

1.  Introduction  and  Background 

[2]  The  flux  of  particulate  and  dissolved  material  from 
New  Guinea  is  very  high  [Mil liman  et  al. ,  1999;  Martin  et 
al.y  2000]  and  it  thus  serves  as  an  excellent  location  to  study 
the  dispersal  of  terrigenous  material  into  the  ocean.  The 
Gulf'  of  Papua  (GoP)  (Figure  1)  is  receiving  a  substantial 
portion  of  this  flux,  as  evidenced  by  its  actively  prograding 
clinoform  [Mil liman  et  al .,  2004].  The  biogeochemical 
processes  that  determine  the  fate  of  terrigenous  sediment 
within  the  GoP  have  been  examined  by  numerous  workers 
[sec  Aller  et  al .,  2004]  as  has  the  distribution  of  trace  metals 
[ Apte  and  Day ,  1998;  Haynes  and  Kwan ,  2002].  These 
studies  indicate  that  terrigenous  material  is  transported  by 
both  coastal  and  deep-water  currents. 

[3]  The  physical  oceanography  of  the  GoP  has  been 
described  in  several  studies  [  Wo  la  ns  hi  et  a  l .,  1984,  1995, 
1999].  The  Coral  Sea  Coastal  Current  (CSCC)  forms  a 
clockwise  gyre  that  is  strongest  during  the  trade-wind 
season.  Across-isobath  tidal  currents  reach  1  m  s_1  in 
response  to  water  level  fluctuations  of  1  -5  m  at  the  mouth 
of  the  Fly  River.  Sub-tidal  along-shclf  currents  of  0.1  m  s^1 
are  coherent  with  the  wind  at  periods  of  about  10  days  on 
the  inner  shelf,  whereas  low-frequency  currents  within  the 
mixed-layer  are  generally  to  the  southwest.  Some  Fly  River 

'Naval  Research  Laboratory,  Stennis  Space  Center,  Mississippi,  USA. 

"Department  of  Geosciences,  Pennsylvania  State  University,  University 
Park,  Pennsylvania,  USA. 

Copyright  2006  by  the  American  Geophysical  Union. 
0094-8276/06/2005GL0254 1 6$05.00 


plume  water  may  infrequently  enter  the  Torres  Straits 
because  of  sub-tidal  water  elevation  differences  between 
the  GoP  and  the  Arafura  Sea.  Modeling  studies  have  been 
used  in  conjunction  with  measurement  programs  to  eluci¬ 
date  dispersal  mechanisms  for  noncohesivc  sediments  in  the 
GoP  [  Wo  l a  ns  ki  et  al .,  1984,  1995;  Hemer  et  al .,  2004]. 

[4]  This  study  complements  earlier  work  by  examining 
the  dispersal  of  dissolved  and  fine-grained  suspended  ter¬ 
rigenous  material  within  the  GoP  for  a  year-long  interval 
using  a  combination  of  hydrodynamic  and  Lagrangian 
tracer  modeling.  Furthermore,  because  of  the  large  model 
domain  used  in  this  study,  we  arc  able  to  study  the  impact  of 
large-scale  flows  on  transport  patterns  without  imposing 
artificial  boundary  conditions.  The  use  of  tracers  to  examine 
complex  multidisciplinary  problems  in  oceanography  has 
led  to  understanding  the  behavior  of  both  suspended  and 
dissolved  constituents  in  a  range  of  ocean  environments. 
[Nakano  and  P ovine c\  2003;  McManus  and  Prandle ,  1997; 
Tomczak  and  Herzfeld,  1998].  Remobilization  processes 
such  as  resuspension  by  waves  have  been  examined  by 
Hemer  et  al.  [2004]  and  are  not  included  in  this  study. 


2.  Hydrodynamic  Simulations 

2.1.  The  EAS16  Model 

[5]  The  Navy  Coastal  Ocean  Model  (NCOM)  is  a  three- 
dimensional  primitive  equation  hydrodynamic  model  [Morey 
et  al .,  2003].  Regional  circulation  during  2003  is  simulated 
using  the  East  Asian  Seas  implementation  of  NCOM 
(EAS16)  with  a  horizontal  cell  size  of  1/16°  (Figure  1 ).  The 
model  assimilates  monthly  mean  discharges  for  51  rivers, 
synthetic  temperature  and  salinity  profiles,  and  altimetry. 
Surface  boundary  conditions  consist  of  heat  fluxes,  precipi¬ 
tation,  surface  pressure,  and  wind  stresses  from  an  operational 
weather  model  (Navy  Operational  Global  Atmospheric  Pre¬ 
diction  System).  The  open  boundary  conditions  include  sea 
surface  elevation,  transport,  temperature,  salinity  and  currents 
provided  by  the  global  NCOM  model. 

2.2.  Circulation  Within  the  Gulf  of  Papua 

[6]  The  typical  surface  circulation  from  the  EAS16 
model  during  the  northwest  monsoon  (Figure  2a)  reveals 
the  clockwise  flow  of  the  CSCC  and  a  weak  southward  flow 
(~0. 1  m  s-1)  on  the  northwest  shelf  as  reported  by  Wolanski 
et  al.  [1995].  The  currents  at  the  entrance  to  Torres  Strait 
have  not  been  reported  for  the  monsoon  but  the  predicted 
surface  currents  of  approximately  0.4  m  s-1  arc  probably 
high,  based  on  observed  oceanographic  features  [  Wolanski 
et  al .,  1984].  The  relative  magnitudes  and  directions,  which 
are  correct,  reveal  highly  variable  flow  due  to  uncorrelated 
sea  level  fluctuations  within  the  GoP  and  Arafura  Sea. 
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Figure  1.  (a)  Bathymetry  used  for  the  EAS16  model.  The 
grid  is  approximately  1/16°  with  a  horizontal  resolution  of 
6-9  km.  The  box  indicates  the  location  of  Figure  lb. 
(b)  Detailed  view  of  the  Gulf  of  Papua  and  northern  Coral 
Sea.  Key:  FR  =  Fly  River;  PR  =  Purari  River;  TS  =  Torres 
Straits;  AS  =  Arafura  Sea;  GP  =  Gulf  of  Papua.  The  stars 
indicate  where  tracers  are  released. 

Predicted  surface  currents  of  0.1 -0.2  m  s-1  are  eastward 
along  the  northeast  coast  of  the  GoP  while  bottom  currents 
(Figure  2b)  are  on  the  order  of  0.01  m  s-1.  Bottom  flow 
varies  in  direction  because  of  the  variable  winds  during  the 
northwest  monsoon. 

[7]  During  the  southeast  trade-wind  season,  the  surface 
currents  on  the  shelf  near  the  Fly  River  vary  from  west  to 
southwest  with  speeds  of  0. 1  0.2  m  s  1  (Figure  3a),  whereas 
bottom  currents  of  less  than  0.1  m  s  1  are  predominantly 


southwest  (Figure  3b)  but  with  more  variability  than  at  the 
surface.  The  surface  (low  along  the  northern  coast  is  weak  and 
variable.  Uniformly  southeast  surface  currents  near  the  north¬ 
east  coast  reach  1  m  s  1  and  the  bottom  llow  is  nearly  aligned 
with  surface  flow.  These  results  are  consistent  with  observa¬ 
tions  [ Wolanski  et  al. ,  1 995 j.  The  modeled  CSCC  forms  a 
closed  clockwise  loop  during  the  trade-wind  season  but  it  has 
less  impact  on  bottom  flow  at  the  shelf  break  than  during  the 
monsoon.  The  predicted  flow  into  Torres  Strait  is  higher  than 
observed  [Wolanski  et  al .,  1 999]  because  the  frictional  effects 
of  several  large  reefs  within  the  northern  Great  Barrier  Reef 
arc  absent  in  the  EAS16  model. 

3.  Tracer  Simulations 

3.1.  The  Particle  Tracking  Model 

[s]  The  particle  tracking  model  uses  three-dimensional 
current  fields  (u.v.w)  and  Gaussian  random-walk  dispersion 
[Lu  et  al .,  2002].  A  particle’s  location  (jr.v.z)  after  a  time 
step  A /  is  given  by: 


-*/+i  —  s,  T  (u  •  A/)  -F  (A Du  •  u/s)  +  (A Du  •  v/s)  +  A Dho  ( 1 ) 

>7+1  —  >7  +  (v  A /)  +  (AD/j-  •  v/s)  +  (A Dlt  •  u/s)  +  A />//0  (2) 
-/+ 1  —  -/  4*  (w  •  A/)  +  ADj  o  (3) 


A 

Z  =  0.0000000:  20030201 0000 
Vector  Scale  =  40.0  cm/s 
140'  144*  148’ 


140'  144'  148' 

B 

Z  =  Bottom:  200302010000 
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Figure  2.  Currents  predicted  by  the  EAS16  model  during 
the  northwest  monsoon  period  in  the  Gulf  of  Papua,  (a) 
Surface  currents,  (b)  Bottom  currents.  Depths  below  100m 
are  shaded. 


Figure  3.  Currents  predicted  by  the  EAS16  model  during 
the  southeast  trade  wind  period  in  the  Gulf  of  Papua,  (a) 
Surface  currents,  (b)  Bottom  currents.  Depths  below  100  m 
are  shaded. 
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where:  it  =  .v-directed  velocity;  A DLT  horizontal  disper¬ 
sion;  s  current  speed;  v  =  v-directed  velocity;  \DHl)  = 
neutral  horizontal  dispersion;  A Dy0  =  neutral  vertical 
dispersion;  iv  =  vertical  velocity.  The  dispersion  terms  are 
given  by: 

A D/t  =  (6  •  Kd i  •  Ar)0'5  •  RAND'  ,  (4) 

AD//o  -  (6  •  Km  •  A/)0.5  •  RAND1  ,  (5) 


AD,  o  =  (6  •  Km  A/)0.5  •  RAND1  ,  (6) 


where  KD],  KD2 ,  and  /fDrare  dispersion  coefficients  with 
values  of  0.002,  0.002,  and  0.0002  m2/s,  respectively.  The 
dispersion  terms  are  randomized  for  fractional  values  from 

1  to  I  using  the  RAND1  j  term,  which  has  a  Gaussian 
distribution.  Particles  arc  introduced  at  specified  points 
using  a  fixed  production  rate  of  2  particles  per  day.  As  a 
particle  settles  through  the  water  column,  equations  ( 1 )  (6) 
are  solved  at  its  new  level  before  the  settling  velocity  is 
applied.  When  the  particle  reaches  the  bottom  it  is  deposited 
and  can  not  be  resuspended. 

[9]  Terrigenous  input  to  the  GoP  is  simulated  using 
tracers  with  settling  velocities  of  10  5  m  s_1  and  10  6  m 
s  '.  The  faster-settling  particles  used  in  this  study  are  more 
representative  of  natural  flocculated  particles  [Fugate  and 
Friedrichs ,  20021  whereas  the  slower  settling  particles  are 
included  to  represent  unflocculated  clay  particles,  woody 
debris,  plankton,  and  dissolved  substances. 

3.2.  Potential  Transport  Paths 

[10]  The  tracer  particles  are  released  at  two  locations 
along  the  northwest  coast  of  the  GoP  (see  Figure  lb).  The 
southern  release  point  is  at  the  mouth  of  the  Fly  River  and 
the  northern  is  west  of  the  Purari  River  in  an  area  with 
several  sources  of  fresh  water.  These  will  be  referred  to  as 
Purari  particles.  The  heavier  particles  (indicated  by  yellow 
circles  in  Figure  4)  settle  near  the  sources  but  the  Purari 
particles  are  more  widely  dispersed  because  of  the  stronger 
and  more  uniform  currents  along  the  northwest  coast.  The 
coastal  currents  transport  the  lighter  particles  eastward,  such 
that  Fly  and  Purari  particles  become  mixed  along  the 
northern  coast.  Some  of  the  Purari  particles  arc  transported 
as  far  as  149°  East.  These  transport  pathways  persist  until 
mid-March. 

[11]  The  transition  to  the  trade  wind  season  occurs  in 
March  and  April,  by  which  time  the  lighter  particles  from 
the  Fly  River  have  reached  the  Purari  but  with  an  offshore 
displacement.  By  mid-March  the  coastal  currents  become 
westward  at  both  rivers  in  response  to  early  intensification 
of  the  southeast  trade  winds.  The  potential  transport  paths 
develop  a  stronger  offshore  component  and  the  lighter 
particles  reach  the  20  m  isobath.  Eastward  transport  along 
the  coast  ceases  by  April  but  some  of  the  lighter  particles 
are  entrained  by  eddies  and  transported  beyond  the  shelf 
break.  Some  of  these  particles  arc  then  transported  eastward 
in  the  CSCC.  The  weakening  of  the  eastward  coastal  current 
is  accompanied  by  an  increase  in  deposition  of  lighter 
particles  along  the  northeast  coastline. 


141°  145°  149° 


141°  145°  149° 


Figure  4.  The  distribution  of  particles  after  one  year  of 
integration.  Topography  is  overlaid  on  model  results  using 
MODIS  imagery.  The  yellow  particles  represent  clay 
minerals  and  the  blue  represent  finer  particles. 

[12]  The  trade  wind  circulation  pattern  becomes  well 
established  by  May  and  lighter  particles  from  the  Fly  River 
enter  the  Torres  Strait  along  its  northern  boundary.  Note, 
however,  that  none  of  the  heavier  particles  enter  the  strait. 
As  the  trade-wind  season  continues,  transport  of  lighter 
particles  through  the  strait  becomes  more  episodic  in 
response  to  sea  level  variations  within  the  GoP  and  the 
Arafura  Sea.  Eastward  transport  of  Purari  particles  contin¬ 
ues  because  of  the  persistent  CSCC  flow  and  these  particles 
arc  entrained  in  eddies  within  the  northern  Coral  Sea.  The 
heavier  particles  from  both  sources  are  deposited  very  near 
the  coast  because  of  the  westward  and  southwestward  shelf 
flow. 

4.  Discussion  and  Conclusions 

[13]  The  final  particle  distribution  reveals  transport  paths 
for  the  faster  settling  particles,  representing  clay  minerals, 
that  are  confined  to  the  inner  shelf  (yellow  circles  in 
Figure  4).  These  potential  transport  paths  are  in  agreement 
with  field  studies,  which  show  high  sedimentation  rates  on 
the  inner  shelf  [Walsh  et  al. ,  2004 1  and  very  little  sediment 
transport  into  Torres  Strait  [Wolanski  et  al ,  1999].  Thus,  it 
appears  that  clays  are  stored  on  the  inner  shelf  along  with 
the  coarser  silts  and  sands  before  being  resuspended  and 
transported  to  the  clinoform  face.  Eastward  transport  of 
clays  could  have  contributed  to  development  of  a  clinoform 
within  the  northeastern  GoP  [Droxler  et  al. ,  2004]. 

[14]  The  lighter  particles  represent  river  input  of  un floc¬ 
culated  clay  minerals,  CDOM,  dissolved  metals,  and  woody 
material.  These  particles  are  widely  dispersed  (blue  circles 
in  Figure  4).  Their  distribution  is  consistent  with  transport 
of  lobster  larvae  by  eddies  within  the  northern  Coral  Sea 
[Dennis  et  al. ,  2001].  Their  transport  into  the  Torres  Straits 
is  in  agreement  with  measurements  of  dissolved  metals 
[Haynes  and  Kwan ,  2002;  Apte  and  Day ,  1998].  The 
potential  transport  pathways  of  dissolved  material  into  the 
GoP  and  Solomon  Sea  to  the  east  are  also  consistent  with 
observations  [Mi l liman  et  al. ,  1999;  Martin  et  al .,  2000]. 
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